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NON-HOMOGENEOUS BOUNDARY VALUE PROBLEMS FOR 
FRACTIONAL DIFFUSION EQUATIONS IN L 2 -SETTING 

KENICHI FUJISHIRO 


Abstract. In the present article, we study the diffusion equations with fractional time 
derivatives. The aim of this paper is to investigate the best possible regularity for the initial 
value/boundary value problems with non-homogeneous Dirichlet boundary data. The main 
tool we use here is called the transposition method. 


1. Introduction 

Let be a bounded domain of M, d with C 2 boundary T := dfl and set Q := Q x (0, T) 
and E := T x (0, T). We consider the following initial value/boundary value problem for a 
partial differential equation with the fractional derivative in time t: 


<9“w + Au = 0 

in Q, 


u = g 

on E, 

(1.1) 

o 

o 

in fl 



with 0 < a < 1. Here df denotes the Caputo derivative, which is defined by 

1 /'' f)ii 

d?u{x,t) = — --/ (■ t-T)- a —(x,T)dr , 

r(i-a)y 0 dr 

and T(-) is the Gamma function (see Podlubny [4]). The differential operator A is given by 
Au(x) = — + c(x)u(x), 

i,j =1 * 1 ' 3 / 

and the coefficients satisfy the following: 

d 

a i:j = dji e C 1 ^), l < i,j <d, ^2 > d\t,\ 2 , x e O, (el 11 , 

i,j =1 

c £ C( H), c(x) >0, x £ Q, 

where /i > 0 is constant. The function g is given on E. 

In the present paper, we study the regularity of the solution to (1.1) in the sense of 
Sobolev spaces. As for this problems, Lions and Magenes [3] showed the result for the 
parabolic equations. 
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2. Main Result 

In this section, we prepare the notations and state our main results. We denote by H S (Q), 
s > 0, the Sobolev spaces. For r, s > 0, we abbreviately set 

H r ’ s (Q) := L 2 ( 0, T; H r (Q)) n H s ( 0, T; L 2 (R)). 

Then H r,s (Q) is a Hilbert space with the norm || • ||#r,s(Q) given by 

II ml H r ’ s (Q) := II L 2 (0,T ;H r (Q)) + IMIiT>(0,T;L 2 (fi))' 

In particular, 

H°'°(Q) = L 2 (Q) = L 2 (0,T;L 2 (R)). 

Similarly we set 

H r '*{ E) := L 2 (0, T ; H r {T)) D H s ( 0, T; L 2 (T)), 

H r ; o S (E) := T 2 (0, T; LT r (r)) D R 0 S (0, T; L 2 (T)). 

with the norm defined by 

II II 2 — II II 2 I II || 2 

Imllimys) •“ ll u llL 2 (o,T;H r (r)) + ll w ll.ff s (o,T;L 2 (r))- 

We clearly have 

H°’°( E) = L 2 (E) = L 2 (0, T ; L 2 (T)). 

As for the spaces with negaive exponents, we dehne 

:= (/rys))', r,s>0. 

The duality paring between R~ r, ~ s (E) and H r ’^(T,) is denoted by {ip,u) rs for x/j G R” r ’ _s (E) 
and u G H r ’^(T). 

We dehne the operator d VA : H s (ti) —> H S ~ 3 / 2 (T), s > 3/2, as 

where i/(x) = (i/i(x),..., Vd{ x )) is the outward unit normal vector to T at x. Then for 
u G H r ' s (Q) with r > 3/2, the trace theorem (Theorem 2.1 in Chapter 4 of Lions and 
Magenes [3]) yields 

#L 6 if-(£), E = 'L = ’U01 ( „ = 0 if s = o) ( 2 . 1 ) 

az /4 r s r 


-p— < C\\u\\ H r,s(Q). 

UL/ A h^CS) 

In order to dehne the weak solution of (1.1), we introduce the dual sytem ; 

{ Hfu + Av — f in Q, 
v — 0 on E, 

4l a u(-,T) = 0 in R, 
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where Df is the backward Riemann-Liouville fractional derivative, which is dehned by 

1 d rT 


D?h(t ) = 


T(1 — a) dt 


(r — t ) a h(r)dr, 0 < a < 1. 


Moreover I!f_ denotes the backward integral of order z/; 


Itf_h{t) = 


(r — t) v 1 h(r)dr, v > 0. 




By the same argument as Chapter 4 in Bajlekova [1], for any / G L 2 {Q ) there exists a unique 
solution v G H 2 ’ a (Q ) of (2.3) such that 


MI// 2 ’“(Q) ^ C1I/IU 2 (<2)- 


(2.4) 


Henceforth we will denote this solution by Vf. Now we apply (2.1) to Vf G H 2,a (Q) and 
obtain 


dv 


dv A 

Since 0 < a < 1, we have 

(see (B.l)). That is, 


L £ H l/2,a/4^ and 


dv 


dv A 


l /!/2,a/ 4 (S) 


< ^\\ V f\\H 2 - a (Q)- 


(2.5) 


H^\ 0, T; L 2 (T)) = H a/ \ 0, T; L 2 (T)) 


9vf 

dis A 


G ih 1/2 ’ a/4 (S) = ih V 0 2 ’ a/4 (S). 


Now we are ready to define the weak solution of (1.1); 


Definition 2.1. A function u is a weak solution of (1.1) if 

(u, f) L 2 (Q) + (g, 9 VA v) 1/2 a/A = 0 (2.6) 

holds for any / G L 2 (Q). 


The main result of this paper is as follows; 


Theorem 2.1. Let g G T 2 (£), then (1.1) has a unique weak solution u G i/ 1 / 2 ’"/ 4 ^) 
satisfying 

\\u\\ H \/2, a / A ( Q } < C\\g\\ L 2(Y,y (2.7) 


Now we roughly describe the strategy of the proof. It is not difficult to show the unique 
existence of the weak solution of (1.1) for g G L 2 (£), but the regularity of cannot 

be directly deduced. Therefore we first show the follwoing two results; 

(i) Regularity of the solution for g G 77 _1//2,_a / 4 (H). 

(ii) Regularity of the solution for g G 7/ 3//2,3 "/ 4 (£). 
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After showing (i) and (ii), we obtain the regularity for g G L 2 (E) by interpolating the above 
results. 

The result for (i) can be easily shown. Indeed, from this definition, we can immediately 
deduce the following proposition; 

Proposition 2.2. Let g G , then (1.1) has a unique weak solution u G L 2 (Q) 

satisfying 

IMU 2 (Q) < (2-8) 

Proof. Combining (2.4) and (2.5), we obtain 


Thus the mapping 

l 2 (Q) 3 / ^ fp- e ff 1/2 '" /4 (S) 

av A 

is bounded, and so is 

L 2 (Q) 3 f — (g, d UA Vf) G C. 

Therefore the Riesz’s representation theorem yields the unique existence of u G L 2 (Q) such 
that 

(«> Dlhq) = - (g, d UA Vf) 

holds for any / G L 2 (Q). Thus we have proven the unique existence of weak solution. 
Moreover for any / G L 2 (Q) we have 

\(u, /)l 2 (Q)| = I {g,dv A Vf) I < ||5 , ||/f-l/ 2 .-“/4(E)ll'9^' i; /ll// 1 /2,a/4 (S ) < C||g|| H -l/2,-a/4( E )||/|| L 2(Q), 

where we have used (2.9) in the last inequality. Therefore we have (2.8). □ 

Thus we have proved part (i). In the next section, therefore, we will consider case (ii). 

3. Regular solution 

We first formulate the functions in hP(0,T) vanishing at t — 0. Following (2.10.3.1b) 
in Triebel [6], for M + := (0, oo) we denote by H S (M. + ) the functions in H S (J8L) which are 
identically zero outside M + ; 

H S (R + ) := {u G H S (R ); suppw C M + } . 

Moreover we set 

S„* + (0,T) :={«|(„, t)\ H‘( R + )}. 

Then similarly to (B.6), we can see that such function spaces have a good property for 
interpolation; 

[H s 0 f(0,T),H°l(0,T)] e = <-' )si+fo2 (0,T), 0 < 5l < s 2 < oo, 0 < 0 < 1. 


<^II/IU 2 (Q). (2.9) 

_ffl/ 2 .“/4(S) 
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We also note that we have the representation of Hq + (0, T) as 

Hq + (0,T) = G H s (0,T)-, £\uW*<oo), 0 < s < 1. 

In particular, for 0 < s < 1, s ^ 1/2, we also have 

\H s (0,T), 0<s<~, 

tfo+(0,T)=| ! 2 

^{u G FP(0,T); u(0) = 0}, -<s<l 

(see (B.l) and (B.5)). For simplicity, we set 

Hg(Q) := L 2 (0, T; H r (Q)) fl H° + ( 0, T; L 2 (F!)) 

and 

^o+(S) := L\ 0, T; IT(r)) D H s 0+ ( 0, T; L 2 (r)). 

As for these spaces, we have the following properties for trace; 

Proposition 3.1. Let u G Hq+(Q), 0 < a < 1, then we have 

i/I e rr3/2,3a/4/y,\ 

u Is e ut 0 + (A). 

Moreover the mapping 

Hq+(Q) 3u^u\ E e H 3 0 i 2 ’ 3a/ \ E) 

is a continuous surjection. 

For the proof of this proposition, we prepare the following lemma; 

Lemma 3.2 (Trace Theorem). Let X and Y be Hilbert spaces such that X is 
to Y densely and continuously. If u € L 2 (M + ;X) D H r (R + ;Y) with r > 1/2, then 

u b \ 0) G [Y, X:}i-(j +1/ 2 ) /r , 0 <j<r- 1/2. 

Moreover, the mapping 

L 2 (M + ; X) n H r (R +] Y) 3u^ u b \ 0) G [Y, X] 1Hj+l/2)/r 
is a continuous surjection. 


(3.1) 


(3.2) 


embedded 


For this lemma, see Theorem 4.2 in Chapter 1 of [3]. 
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Proof of Proposition 3.1. Without loss of generality, we may consider the case of 

Q = = {( Xl ,..., x d -i,x d y, xi,..., x d -i £l, x d > 0}, 

r = M 0 '- 1 = {(xi, ..., x d - 1 ); xi, ..., x d -i G R}. 

We abbreviately write x' := (xi,... ,x d _ i). Then we have 

U 6 = £ 2 (0, T; 2 (Rp) n HJ + (0, T; L 2 (Rp) 

if and only if 

u 6 L 2 x (0,X))) n H 2 (R +A ,; L 2 (Rjr‘ x (0,T))) 

= I 2 (R+, ff 0 2 i”(S)) n H\ R + , L 2 (E)). 

We apply Lemma 3.2 as X — H^(E) and Y — L 2 (£). Then we have 

«b = «(0) € [L 2 (S),2/ 2 + “(E)] 3/4 = H 2 i 2M/i ( E). 

The surjectivity also follows from Lemma 3.2. □ 


By using the trace theorem stated above, problem (1.1) with g G Hy+' :in - 4 (Y J ) is directly 
reduced to the following problem with homogeneous boundary condition; 


{ d^u + Au = F in Q, 

u = 0 on E, (3.3) 

u(-, 0) = 0 in fl. 


For (3.3), Gorenflo, Luchko and Yamamoto [2] showed the L 2 -maximal regularity. In their 
setting, the Caputo derivative 9“ equipped with the initial value 'u(O) =0 is formulated as 
an operator in L 2 (0, T ) with its domain given by 


V{d?) = 


H a (0, T), 
u e !L 1 / 2 (0,T); 


,2 dt 


I \u{t)Y — < oo 
0 t 

{u G H a ( 0, T); w(0) = 0} , 


0 < a < 1/2, 
a = 1/2, 

1/2 <a< 1 . 


(3.4) 


By (3.1) and (3.2), this can be rewritten as 

D(a?) = S 0 « + (0,T). (3.5) 

Thus we can see that if (3.3) has a “solution” in T?(<9"), then the initial condition u(-, 0) = 0 
is satisfied in a weaker sense. They also revealed that the above operator is essentially 
equivalent to the Riemann-Liouville derivatives, which were already discussed in [1]. Anyway 
we obtain the following result; 


Lemma 3.3. Let 0 < a < 1 and F G L 2 {Q), then (3.3) has a unique solution u G Hq’/*(Q) 
satisfying 

IMI H^{Q) < C\\F\\ l i {Q) . 
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For the proof of this lemma, see Theorem 4.3 in [2], 

Proposition 3.4. Let 0 < a < 1 and g £ //q| 2 ’ 3q ^ 4 (S). Then problem (1.1) has a unique 
solution u £ H 2 f(Q) satisfying 

IMI T7 2 '“(Q) < C , ||fl'||tf3/2,3 Q ,/4( E ). 


Proof. Since g £ H ( 


3/2,Za/A 

0 + 


(£), Proposition 3.1 yields that there exists g £ Hq+(Q) such that 


By setting F : = 
w £ Hlff(Q) of 


which satisfies 


9 Is — g and \\g\\ H 2 , a ^ < C\\g\\ H 3/2, 3 a/A^y 

—Ag — dfg £ L 2 (Q ) and applying Lemma 3.3, there exists a solution 


dfw + Aw = F 

in Q, 

w = 0 

on £, 

o 

o 

in fl, 


IMI H 2 ’“(Q) < C\\F\\ L 2 {Q) < C\\g\\ H 2 , a (Q) < C\\g\\ H 3/2,3 a /4^y 
Then u := w + g satisfies (1.1) and 

IMI H 2 ’ a (Q) < IMM“(Q) + IMM 2 ’“(Q) < C , ||5 , ||7/3/2,3a/4( S ). 

Thus we have completed the proof. □ 


4. Proof of the main result 

In this section, we complete the proof of Theorem 2.1 by interpolation. 


Proof of Theorem 2.1. Let n be the operator which operates the boundary data g to the 
weak solution u of (1.1). Then, by Propositions 2.2 and 3.4, we have 

7T e £ ( H -VW 4 ( S ) ; L 2 (Q)) n C (tf 3 { 2 ’ 3 “ /4 (£); tf 2 ’“(Q)) , 

where £(X, Y ) denotes the set of linear and bounded operators from X to Y. By Proposition 
A.l, the operator tt also belongs to 

C ([tf-‘/ 2 '-“/ 4 (E), E)]»; [L 2 (Q), iA“(Q)]„) 

for any 0 < 6 < 1. In particular, if we set 6 = 1/4, then 

M”i/2’-“/4(S), tf 3 { 2 ’ 3a/4 (£)] 1/4 = L 2 (£) and [L 2 (Q), H 2 ’ a (Q)} 1/4 = E x ' 2 ^\Q) 1 

and therefore we have 

7T £ £(L 2 (S);LT 1/2 ' Q/4 (g)). 

Thus we have completed the proof. □ 
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Appendix A. Interpolation 

Throughout this article, we often use the word “interpolation” as a complex interpolation 
defined bellow. As for the detailed argument on this topic, we can refer to Triebel [6], 
Yagi [7] and the references therein. On the other hand, in some classical works such as 
Lions and Magenes [3], the “interpolation” of two Hilbert spaces is defined as the domain of 
fractional powers of positive and self-adjoint operator. We will see that these two kinds of 
definitions coincide with each other (see Proposition A.2). Therefore, we can refer to [3] and 
use some of their results (e.g., Theorem 4.2 in Chapter 1 of [3]) without any confusion. In this 
section, we recall the definition of complex interpolation of Banach spaces and summarize 
their fundamental properties. 

Let X t be a Banach space equipped with the norm || • || Xt (* = 0, 1) and suppose that Ad 
is embedded in A 0 continuously and densely. Let S be defined by 

S := {z E C; 0 < Rex; < 1}. 

We say that a function F : S —> X 0 belongs to H(X 0 ,Ad) if and only if the following 
conditions (H1)-(H3) are satisfied; 

(HI) F is analytic in S. 

(H2) F is bounded and continuous in S. 

(H3) R 3 y i —y F( 1 + iy) E Ad is bounded and continuous. 

It is known that TL(X 0 , X]) is a Banach space with the norm || • ||% given by 

\\F\\n ■= max ("sup \\F(iy)\\ Xo , sup \\F(1 + ij/)||x 1 ') , F E H{ Ad,Ad). 

\yeR yeR J 

For each 0 < 6 < 1, we define the space [A 0 , Ad]y by 

[A 0 , Xi]g := {u E X 0 ] u = F(6) for some F E H{X 0l Ad)} 

Moreover [Ao, Ad }q is a Banach space with the norm || • ||# defined by 

\\u\\ e := inf ||F|| W , uElX^X^. 

FeH(X 0 ,X i), 

F(8)=u 

By the interpolation, we can show various kinds of “intermediate properties”. For example, 
if a linear operator T is bounded from A" 0 into Yo and from Ad into Y\ at the same time, 
then we can deduce that T is also a bounded operator from [A 0 , X\\o into [Yj, Y]o for any 
0 < 6 < 1 . 

Proposition A.l. Let Ad (resp. Y 1 ) be embedded to X 0 (resp. Y 0 ) densely and continuously. 
Then for any 0 < 6 < 1, 

C(X Q , Y 0 ) C £(Ad, Yi) C C{[X 0, Ad] 0 , [Y 0 , Y 1 ] e ) 

and we have 

l|r|L(|x 0 .x I ]»,|v 0 ,y ll .) < T e C( A'o, Y„) n £(X„n). 
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Moreover we can also characterize the domain of fractional power of operators; 


Proposition A.2. Let X be a Hilbert space and A : X —y X be a positive and self-adjoint 
operator. Then we have 

V(A e ) = [X,V(A)} e , 0<6<1 

with isometry. 


Here we note that [A, U(A)]g stated above coincides with [D(A), X]i_ e in the notation by 
Lions and Magenes [3]. 


Appendix B. Sobolev spaces 


Let be a domain of W l with smooth boundary T = dfl. We dentoe by H S (Q), s > 0, the 
space of Bessel potentials (see e.g., [6]). We can see that H s (Ll), s > 0, has a good property 
with respect to the interpolation; 

[H Sl (n),H S2 (Q)}e = H^~ e)si+dS2 (n), 0 < Sl < s 2 < oo. 


As a characterization of subspaces of LP(fl) consisting of the functions vanishing on T, we 

often use Hq(Q) (also denoted by H s (fl)), the closure of C^°(r2) in H S (Q). We note that 
Hq(Q) for 0 < s < 1 has the representation as 




1 [u G H s ( fl); n|r = 0}, - < s < 1. 


(B.l) 


However, we see that for these spaces, good properties concerning the interpolation fails 
for some “singular” cases—when s = integer+1/2. Indeed, according to Theorem 11.6 in 
Chapter 1 of [3], the identity 


is valid for the case in which si,S 2 , (1 — 0)si + ds 2 ^integer+1/2. On the other hand, by 
Theorem 11.7 in Chapter 1 of [3], if (1 — 9)si + 0s2 = n + 1/2 for some integer // (si, S 2 are 
still assumed not to be integer+1/2), then 

K +1,2 ( Si) := (B.2) 


In particular, for fi — 0 we have 

[z. 2 (n),//i(n)] 1/2 = i/+(n) £ h'J 2 ( si) = /A 2 (Si). 

We note that the space also can be characterized without interpolation. In fact, 

by (11.52) in Chapter 1 of [3], we have the repredentation as 

H^ l/2 {Tl) = |ue Hi; +1/2 (Q)- p- 1/2 dfu G L 2 (H) for |a| = /x} , (B.3) 
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where p : 0 —* [0, oo) is a smooth function such that 


lim 


p(x ) 


= d ^ 0 , Jo 6 r. 


x^-xo d^x, r) 

Moreover by substituting p — 0 in (B.3), we can rewrite (B.3) as 

\u(x)\ 2 


4 / 2 ( 0 )= 




P(x) 


-dx < oo 


Thus we see that each element in H^ 2 ^) vanishes on T in a certain sense. 

In the following, therefore, we introduce another formulation of functions in H S (Q) vanish¬ 
ing on T, which includes Hqq 1 ^ 2 ^) as a particular case. We denote by H S (Q) the subspace 
of H s (R d ) consisting of the functions which are identically zero outside O (see (4.3.2.1b) in 
[6]). That is, 

H s (Q) := {u G H s (R d ); supp u C 0} . (B.4) 

By regarding elements of H S (Q) as functions defined on hi, we have 

H S (Q) = Hq(Q), s > 0, s / integer + -. (B.5) 

Moreover, we have the following good property for interpolation (see Corollary 1.6 in Chapter 
3 of Strichartz [5]); 

[# sl (^)>^ S2 (^)]e = TT (1 - 0)si+0S2 (f2), 0 < si < s 2 < oo, 0 < 6 < 1. (B.6) 


By comparing (B.2) with (B.6) and noting (B.5), we have 

This identity also can be verified by comparing the representations (B.3) and (2.4.2.7) in [6] 
with p = 2. Thus the subspace 7P(fl) introduced in (B.4) is more appropriate than 
when we deal with interpolation. 
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